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Abstract 

A model for total cross sections with virtual photons is presented. In par- 
ticular 7*p and 7*7* cross sections are considered. Our approach extends 
on a model for photoproduction, where the total cross section is subdivided 
into three distinct event classes: direct, VMD and anomalous. With in- 
creasing photon virtuality, the latter two decrease in importance. Instead 
Deep Inelastic Scattering dominates, with the direct class being the 0{a s ) 
correction thereof. Hence, the model provides a smooth transition between 
the two regions. By the breakdown into different event classes, one may 
aim for a complete picture of all event properties. 
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1 Introduction 



Traditionally, different descriptions are used for virtual and real photons. Virtual photons 
in the DIS (Deeply Inelastic Scattering) region are normally described as devoid of any 
structure, while for the real ones, the possibility of hadronic-like fluctuations play an 
important role. In the region of intermediate Q 2 , it should be possible to find a description 
starting from either extreme. Then the language may not always be unique, i.e. a given 
Feynman diagram may be classified in different ways. 

Several paths are possible; one is to have an explicit description in terms of higher- 
order Feynman diagrams 0, which then may contain several hard scales. However, in 
this article a main point is to obtain a smooth transition to soft physics, beyond the 
region of validity of perturbation theory, and then it is convenient to be able to consider 
transitions in one scale at a time. In the following, we will therefore develop one specific 
approach, where the main principle is to characterize events by a set of standard scales, 
such as the photon virtuality (or virtualities for 77), the photon resolution scale(s), and 
the hard scale of a partonic subprocess. Depending on the relative ordering of these 
scales, events are classified in different categories. Special emphasis is normally put on 
the hardest scale of the event. This may determine e.g. whether an event is classified as a 
DIS or a resolved photon one. Matrix elements are defined to lowest required order only, 
and higher-order corrections are approximated by parton showers. For instance, if the 
hard scale is larger than the resolution scale of the photon, a partonic evolution is allowed 
between the two scales. The DGLAP equations are then suitable for this evolution and 
are to be combined with the appropriate leading-order matrix elements, responsible for 
the hard scattering. In a 7p event, the photon resolution scale may well be smaller than 
both the photon virtuality and the hard-scattering scale. In this way, partly unordered 
evolutions in pj_ are accounted for. This is still less complete than allowed by the CCFM 
equations |3], £|], where the ordering may be broken in several steps, but could well provide 
the bulk of non-ordering effects at current energies. 

The classification used here is not an economical route, however, since it leads to many 
event classes. For studies e.g. of the total cross section in the intermediate-Q 2 region, 
it is cumbersome and not necessarily better than existing approaches ||. However, by 
the breakdown into distinct event classes, the road is open to provide a (more or less) 
complete picture of all event properties. It is this latter aspect that has then guided the 
model development. 

As an example, the 7*p cross section is divided into a VMD, an anomalous, a direct 
and a DIS component. In the limit Q 2 — > 0, the DIS process 7*q — > q becomes kinemati- 
cally forbidden and only the first three event classes remain, reproducing the result from 
the photoproduction model (except for some improvements). On the contrary, when Q 2 
increases from zero to high values, the resolved processes decrease in importance (as given 
by dipole dampening factors) , the direct ones also drops (by Q 2 dependence and a shrink- 
ing phase space) and finally only the DIS process remain. At intermediate Q 2 values, the 
direct processes and the DIS (+parton showers) process overlap, since, in some regions 
of phase space, they are equally valid descriptions of the same physics. It thus becomes 
necessary to avoid double-counting, e.g. by introducing Sudakov style form factors for 
the DIS process, suppressing those parton configurations covered by the direct processes. 

Based on the study of jet production by virtual photons || we extend the modeling to 
\ow-p± events. Clearly, a smooth transition from perturbative to non-perturbative physics 
is desired, and is achieved partly by allowing for multiple parton-parton scatterings. These 



1 



Figure 1: Deeply inelastic charged lepton-hadron scattering. (The four-momenta 
are given in parentheses.) 

are needed to explain the underlying event activity seen in the data in hadron-hadron 
collisions J7| , and the ideas will be taken over to the hadronic resolved photon components. 
At the same time, the multiple parton-parton interactions unitarize the jet cross sections, 
a necessity for a decent growth of the total cross section with increasing energy. 

The outline of the paper is as follows. First, the different event classes are presented, 
including an extension to virtual photons of a model for photoproduction. How the 
event classes are combined, in order to avoid double-counting but still to cover the whole 
phase space, is described in detail. Then, results of the model are presented, with some 
comparisons to available data. In particular, total cross sections with virtual photons are 
shown with emphasis on the importance of the different event classes at various photon 
virtualities. Some distributions of event properties are shown to vary significantly with the 
photon virtuality, others less so, illustrating a complex model working smoothly between 
the regions of real and highly virtual photons. Similarities in 7p and 77 events will be 
pointed out where appropriate. Finally, a summary and an outlook is given. 



2 Event Classes 

In the following, we begin by a reminder on the models for DIS and photoproduction, 
before embarking on the generalization also to intermediate virtualities in 7*p processes. 
The 77, 7*7 and 7*7* processes thereafter follow by an application of the same rules. 



2.1 Deeply Inelastic Scattering 

The Deeply Inelastic Scattering of a high-energy charged lepton off a proton target, Fig. [1], 
involves a single electroweak boson exchange between a beam lepton and a target quark. 
At not too large Q 2 only photon exchange need be considered. Then the double- differential 
ep cross-section for DIS can be expressed in terms of the total cross-section for virtual 
transverse (T) and longitudinal (L) photons M: 
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with the fluxes 



/yy.Q 2 ) = (3) 



The conventional DIS variables 



are related e.g. by 

Q 2 = xy2kP, W 2 = (1 -x)y2kP + m 2 p , (5) 

where 2A;P « (/c + P) 2 = s. Thus there are only two kinematical degrees of freedom. 
The cross-sections can be related to the proton structure functions F 2 and F L by |], |TtJ 



r? I r\1\ & X )Q 2 i \ t-, / <n 2\ (1 x)Q 2 

F2i ^ Q } = 47r 2 a cm (Q 2 + 4mgx 2 ) ^ T + ^ ' Q } = 4^(^ + 4^) ^ (6) 

and the total virtual photon-proton cross section by 

7*p_ , 47r 2 a cm (Q 2 + Am 2 x 2 ) 4ir 2 a em 2 

*~ = + = (1 _ g)g2 F 2 (x, Q ) * F 2 (x, Q). (7) 

In the parton model, 

F 2 (x,Q 2 ) = J24 {xq(x,Q 2 )+xq(x,Q 2 )} , F L (x,Q 2 )=0, (8) 
q 

to lowest order. Such an interpretation is not valid in the limit Q 2 — > 0, where gauge 
invariance requires F 2 (x,Q 2 ) —>■ so that remains finite. We will replace the DIS 
description by a photoproduction one in this limit. Hence, at small photon virtualities, 
the DIS process 7*q — > q should be constructed vanishingly small as compared to the 
contribution from the interaction of the hadronic component of the photon (to be discussed 
in the next section). To obtain a well-behaved DIS cross section in this limit, a Q 4 /{Q 2 + 
m 2 ) 2 factor is introduced. Here m p is some non-perturbative hadronic parameter, for 
simplicity identified with the p mass. One of the Q 2 / (Q 2 + m 2 ) factors is required already 
to give finite cr/Jt f° r conventional parton distributions, and could be viewed as a screening 
of the individual partons at small Q 2 . The second factor is chosen to give not only a finite 
but actually a vanishing cr^g for Q 2 — > in order to retain the pure photoproduction 
description there. This latter factor thus is more a matter of convenience, and other 
approaches could have been pursued. Then, in the parton model, eq. (|7|) modifies to a 
DIS cross section 

a ™* wTS^ < M*,q 2 ) + *,q 2 )} • 0) 

For numerical studies, the available parton distribution parameterizations for the proton 
have some lower limit of applicability in both x and Q 2 . For values below these minimal 
ones, the parton distributions are frozen at the lower limits. 
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a) b) c) 

Figure 2: Contributions to hard 7p interactions: a) VMD, b) direct, and c) anoma- 
lous. Only the basic graphs are illustrated; additional partonic activity is allowed 
in all three processes. The presence of spectator jets has been indicated by dashed 
lines, while full lines show partons that (may) give rise to high-p^ jets. 



2.2 Photoproduction 

To first approximation, the photon is a point-like particle. However, quantum mechani- 
cally, it may fluctuate into a (charged) fermion-antifermion pair. The fluctuations 7 <-> qq 
can interact strongly and therefore turn out to be responsible for the major part of the 
7P total cross section. The total rate of qq fluctuations is not perturbatively calcula- 
ble, since low-virtuality fluctuations enter a domain of non-perturbative QCD physics. 
It is therefore customary to split the spectrum of fluctuations into a low-virtuality and 
a high-virtuality part. The former part can be approximated by a sum over low- mass 
vector-meson states, customarily (but not necessarily) restricted to the lowest-lying vec- 
tor multiplet. Phenomenologically, this Vector Meson Dominance (VMD) ansatz turns 
out to be very successful in describing a host of data. The high-virtuality part, on the 
other hand, should be in a perturbatively calculable domain. 
In total, the photon wave function can then be written as 

|7) = Cbare|7bare) + C V \V) + ^ ^1^) + C ^ ^) ( 10 ) 

V=p°,w,(f>,J/ip q=u,d,s,c,b l=e,y,,T 

(neglecting the small contribution from T). In general, the coefficients q depend on the 
scale /i used to probe the photon. Thus cj ~ (a cm /2ii)(2/3) \n(fi 2 /mj). Introducing a cut- 
off parameter ko to separate the low- and high-virtuality parts of the qq fluctuations, one 
similarly obtains Cq ~ (a em /27i)2e^\n(fi 2 /k^). The VMD part corresponds to the range of 
qq fluctuations below k and is thus /^-independent (assuming fi > k ). In conventional 
notation c v = 4na em / fy, with f v /4:7i determined from data to be 2.20 for p°, 23.6 for 00, 
18.4 for (j) and 11.5 for J/ijj |I0[ . The ko parameter is constrained by fits to i 7 ^, i.e. to 
the parton distributions of the photon, to be k ~ 0.6 GeV |TT[. (The fits also contain 
other model uncertainties, and are only logarithmically dependent on k , so the precision 
is not high.) Finally, c bare is given by unitarity: c^ aro = Z 3 = 1 - J2 c v ~ E c l ~ E c v 
In practice, Cb are is always close to unity. Usually the probing scale \i is taken to be the 
transverse momentum of a 2 — > 2 parton-level process. 

The subdivision of the above photon wave function corresponds to the existence of 



three main event classes in 7p events |[L2| , cf. Fig. |2|: 

1. The VMD processes, where the photon turns into a vector meson before the inter- 
action, and therefore all processes allowed in hadronic physics may occur. This in- 
cludes elastic and diffractive scattering as well as low-p_i_ and high-pj_ non- diffract ive 
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events. 

2. The direct processes, where a bare photon interacts with a parton from the proton. 

3. The anomalous processes, where the photon perturbatively branches into a qq pair, 
and one of these (or a daughter parton thereof) interacts with a parton from the 
proton. 

All three processes are of 0(a em ). However, in the direct contribution the photon structure 
function is of 0(1) and the hard scattering matrix elements of 0(a em ), while the opposite 
holds for the VMD and the anomalous processes. The £ + £~ fluctuations are not interesting 
for us, and there is thus no class associated with them. 

The difference between the three classes is reflected in the beam jet structure. The 
incoming proton always gives a beam jet containing the partons of the proton that did not 
interact. On the photon side, the direct processes do not give a beam jet at all, since all 
the energy of the photon is involved in the hard interaction. The VMD ones (leaving aside 
the elastic and diffractive subprocesses for the moment) give a beam remnant just like 
the proton, with a moderately small 'primordial k±' smearing. The anomalous processes 
give a beam remnant produced by the 7 — > qq branching, with a transverse momentum 
going from k upwards. 

Based on the different event classes discussed above, the total photoproduction cross 
section can be written as 

7P _ 7P 1 7P 1 7P CI 1 ) 

°tot — °VMD ^ direct ^ "anomalous ■ 1 11 J 

Total hadronic cross sections show a characteristic fall-off at low energies and a slow 
rise at higher energies. This behaviour can be parameterized by the form 

a AB (s) = X AB s e + Y AB s~ v (12) 



for A + B — > X. The powers e and rj are universal, with fit values [13 



e » 0.0808 , 77 » 0.4525 , (13) 

while the coefficients X AB and Y AB are process-dependent. Equation (|T2| ) can be in- 
terpreted within Regge theory, where the first term corresponds to pomeron exchange 
and gives the asymptotic rise of the cross section. The second term, the reggeon one, 
is mainly of interest at low energies. For the purpose of our study we do not have to 



rely on the Regge interpretation of eq. (12), but can merely consider it as a convenient 
parameterization. 

The VMD part of the 7p cross section is an obvious candidate for a hadronic descrip- 
tion. The diagonal VMD model suggests: 

<md(*)= E ^ (14) 

V=p°,u>,(f>,J/tf? ^ V 

Assuming an additive quark model the V^p cross sections can be parameterized as |12| 

^ot P (s) ~ o%(8) w ~ (ct^p + ct^o?) » 13.63 S e + 31.79s-" [mb] , 
o&(s) w o£t P + o£" p - » 10.01^-1.52*-" [mb], (15) 



m j/V iU 
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with s in GeV 2 . Adding the vector meson contributions, we arrive at 



< MD ( S ) « 53 - 4 * £ + 115s ~" [H • ( 16 ) 

There is no compelling reason that such an ansatz should hold also for the total 7p 
cross section, but empirically a parameterization according to 

a^(s) = X^s e + Y^p s -" = 67.7s e + 129s"" [fib] (17) 



does a good job [|13|]. For instance, these parameter values were used to predict the 
high-energy behaviour of the cross section, close to what was then measured by HI and 
ZEUS. Thus VMD corresponds to approximately 80% of the total 7p cross section at 
high energies, with the remaining 20% then shared among the direct and anomalous 
event classes. 

The anomalous contribution can be written as 

n, r°° Ah 2 

< om ^s) = ^HK L^r *J ( 18 ) 

where the pref actor and integral over dkj_/ k\ corresponds to the probability for the photon 
to split into a qq state of transverse momenta ±kj_. The cross section for this qq pair to 
scatter against the proton, <r qqp , need to be modeled. Based only on geometrical scaling 
arguments (to be discussed in the next section), one could expect a decrease roughly like 
l/kj_. This suggests an ansatz 

k 2 - 

a wp( s; k ± ) = a y(qq)p (s) . (19) 

k{ 

The kv(qq) is a free parameter introduced for dimensional reasons. It could be associated 
with the typical k± inside the vector meson V formed from a qq pair: p° « uj for u and d, 
(f) for s, J/ip for c. As a reasonable ansatz, one could guess fcy( q q) ~ my/2 « m p /2. (For 
heavier quarks, a higher mass scale is indicated, but also a correspondingly larger lower 
integration limit in k^_, so the two effects cancel more or less.) Fits to the total cross section 
at not too high energies, with a large VMD and a small direct contributions subtracted, 
give corresponding numbers, kv(qq) ~ 0.4 GeV for a k ~ 0.5 GeV. These values are 
here related to each other: if the latter were to be changed, the former would have to be 
retuned accordingly. In the following we use this set of numbers. The anomalous cross 
section can thus be written as 

< P o malous (s) - ^ E 2 < |r ^tf ^ )P ^ ( 20 ) 

To leading order, the direct events come in two kinds: QCD Compton 7q — > qg 
(QCDC) and boson-gluon fusion 7g — > qq (BGF). The cross sections are divergent in the 
limit k± — > for the outgoing parton pair. Therefore a lower cut-off is required, but no 
other specific model assumptions. 

2.3 Combining the photoproduction processes 

The VMD, direct and anomalous classes have so far been considered separately. The com- 
plete physics picture presumably would provide smooth transitions between the various 
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Figure 3: (a) Schematic graph for a hard 7p process, illustrating the concept of two 
different scales, (b) The allowed phase space for this process, with one subdivision 
into event classes. 



possibilities. To understand the relation between the processes, consider the simple graph 
of Fig. |3]a. There two transverse momentum scales, k± and p±, are introduced. Here k± 
is related to the 7 — > qq vertex while p± is the hardest QCD 2 — > 2 subprocess of the 
ladder between the photon and the proton. (Further softer partons in the ladders are 
omitted for clarity.) The allowed phase space can then conveniently be represented by a 
two-dimensional plane, Fig. ^]b. The region k± < k corresponds to a small transverse 
momentum at the 7 — > qq vertex, and thus to VMD processes. For k± > ko, the events 
are split along the diagonal k±_ = p±. If kj_ > p±, the hard 2^2 process of Fig. |^a is 
7g — > qq, and the lower part of the graph is part of the leading log QCD evolution of the 
gluon distribution inside the proton. These events are direct ones. If p± > k±, on the 
other hand, the hard process is qq' — > qq', and the 7 — > qq vertex builds up the quark 
distribution inside a photon. These events are thus anomalous ones. 

It should be remembered that the classification is only simple away from the border 
regions. When k± ~ p±, say, a description either in terms of anomalous or direct inter- 
actions would be possible. Also higher-order corrections will blur the picture, although 
pragmatic separations normally can be found. 

A comment on the choice of variables. Instead of the squared transverse momenta 
k\ and p\ one could equally well have imagined the squared virtualities —k 2 and —p 2 . 
However, in many processes several Feynman graphs contribute, both in the t-, u- and 
s-channel. The association of p 2 = t or p 2 = u is then not unique, while p\ = tii/s is 
well-defined and coincides with —t or — u in the respective singular limit. Whether events 
are classified by one kind of scale or another should not be relevant, at least so long as one 
is consistent and stays with a leading-order description. Only for the virtual photon do 
we use the scale Q 2 = —q 2 = —t rather than q 2 ±1 in order to stay in line with conventional 
notation, and because there is no scale choice ambiguity here. 

What complicates the picture is that an event may contain several interactions, once 
one considers an incoming particle as a composite object with several partons that may 
interact, more or less independently of each other, with partons from the other incoming 
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particle. Such a multiple parton-parton interaction scenario is familiar already from 
pp physics 0. Here the jet cross section, above some p_i_ m m scale of the order of 2 GeV, 
increases faster with energy than the total cross section. Above an energy of a few hundred 
GeV the calculated jet cross section is larger than the observed total one. Multiple 
interactions offers a solution to this apparent paradox, by squeezing a larger number of 
jet pairs into the average event, a process called unitarization or eikonalization. The 
perturbative jet cross section can then be preserved, at least down to p± m i n , but in the 
reinterpreted inclusive sense. At the same time, the unitarization plays a crucial role in 
taming the growth of the total cross section. 

The composite nature of hadrons also fills another function: it regularizes the singu- 
larity of perturbative cross sections, such as qg — > qg, in the limit p± — > 0. Perturbative 
calculations assume free colour charges in the initial and final states of the process, while 
the confinement in hadrons introduces some typical colour neutralization distance. It is 
the inverse of this scale that appears as some effective cutoff scale p± m i n — 2 GeV, most 



likely with a slow energy dependence [14]. One possible parameterization is 



P±mUs) = (1.9 GeV) (p^f , (21) 

with the same e as in eq. (0), since the rise of the total cross section with energy via 
Regge theory is related to the small-x behaviour of parton distributions and thus to the 
density of partons. 

Now, if an event contains interactions at several different p± scales, standard practice 
is to classify this event by its hardest interaction. Several reasons can be given; one is 
that a softer interaction may be confused with QCD radiation emitted from the harder 
one, and thus cannot be identified on an exclusive basis. With this prescription, the cross 
section for an event of scale p± is the naive jet cross section at this pj_ scale times the 
probability that the event contained no interaction at a scale above p±. The latter defines 
a form factor, related to probability conservation, analogously to but not equivalent with 
the Sudakov form factor of parton showers. At large p± values the probability of having 
an even larger p± is small, i.e. the form factor is close to unity, and the perturbative 
cross section is directly preserved in the event rate. At lower p± values, the likelihood of 
a larger p± is increased, i.e. the form factor becomes smaller than unity, and the rate of 
events classified by this p± scale falls below the perturbative answer. 

We expect this picture to hold also for the VMD part of the photon, since this is clearly 
in the domain of hadronic physics. Thus, in the VMD domain k± < ko, the region of large 
pj_ in Fig. |b is populated according to perturbation theory, though with nonperturbative 
input to the parton distributions. The region of smaller p± is suppressed, since the form 
factor here drops significantly below unity. 

As one moves away from the 'pure' VMD states, such as the p°, much of the same 
picture could well hold. Interactions at a larger k± value could be described in terms 
of some p' state. The uncertainty relation gives us that a state of virtuality ~ k± has 
a maximal size ~ l/k± and thus spans an area oc l/k\. In a naive picture, the p'p 
cross section would then drop with increasing k±, but flatten out at a value given by the 
proton size. We know that this is not the scaling observed in nature, however, where 



e.g. the J/Vp cross section is much below the 0p one |15| . Rather it appears that the 
cross section is proportional to the area of the state interacting with the proton, i.e. a 
(kind of) geometrical scaling. Such a behaviour becomes understandable (but not easily 
predicted in detail) when one remembers that the colour neutralization distance inside a 
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Figure 4: An alternative classification of the phase space in Fig. 
takes into account unitarization effects. 



which better 



more virtual photon state is also reduced, so that the interactions in general tend to be 
weakened by interference effects not included in the simple perturbative cross sections. 
This could then be the origin for a geometrical scaling like the one in eq. ([19]). 

Alternatively to the geometric approach, the anomalous and direct cross sections for 
such a k± state could be calculated. If these latter are larger than the geometrical scaling 
answer, one could expect multiple interaction effects to become important and couple 
the two event classes. That is, the presence of an anomalous interaction of some large 
p± would preempt the event from being classifiable as a direct one with a lower p±, in 
analogy with the discussion on hadronic physics above. One can even attempt eikonalized 
descriptions Hl~6| , although such require a number of assumptions to be made. 

Calculating the perturbative anomalous cross section in the region given by p± > 
max(A;_|_,p^ min (s)), the geometric scaling answer is exceeded for some region kj_ r^i k\^ 
with ki « 2 — 4 GeV (higher for higher energies). Only for k±_ > k\ is the jet cross 
section dropping below the geometric scaling one. At these larger k± values, the direct 
rate dominates over the anomalous. As a convenient but rather arbitrary choice, for 
subsequent studies we put k\ = p± m i n (s), with the latter given by eq. (^1|). This value lies 
below the crossover point noted above, but the simple modified geometrical scaling picture 
is not precise enough that it can be trusted too far, so a separate k\ parameterization 
would seem excessive. 

The final scenario is illustrated in Fig. The bulk of the cross section, in the region 
k± < k\, is now described by the photon interacting as dense, hadronic states, VMD for 
k± ~ k and Generalized VMD (GVMD) for fc ~ k± £ k\. The total VMD cross section 
is given by the pomeron-type ansatz, while the jet cross section can be obtained from the 
parton distributions of the respective vector meson state. Correspondingly, the GVMD 
states have a total cross section based on Pomeron considerations and a jet cross section 
now based on the anomalous part of the parton distributions of the photon. In principle, 
an eikonalization should be performed for each GVMD state separately, but in practice 
that would be overkill. Instead the whole region is represented by one single state per 
quark flavour, with a jet production given by the full anomalous part of the photon 
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Figure 5: (a) Schematic graph for a hard 7*p process, illustrating the concept of 
three different scales, (b) Event classification in the large-Q 2 limit. 

distributions. By extending this jet rate also to the region k± > k%, one does introduce 
a small mismatch between the low-p^ and high-p^ descriptions, but retains the correct 
full jet rate at high p±'s. Implicitly, the presence of a spectrum of k± states can be taken 
into account by the choice of a realistic p± spectrum of the photon beam remnant, maybe 
still without the full correlation structure between the high-p^ and low-p^ parts of the 
individual events. 

Thus, post facto, the approximate validity of a Regge theory ansatz for a^t is making 
sense. Above k% only the direct cross section need be considered, since here the anomalous 
one is negligibly small, at least in terms of total cross sections. (As noted above, we have 
actually chosen to lump it with the other GVMD contributions, so as not to lose the jet 
rate itself.) Numerically, the recipe of extending the anomalous contribution to infinity 
according to a GVMD scaling recipe, as is done in eq. (p0|), is about equally good. The 
latter may involve some double-counting with the direct cross section, but not more than 
falls within the general uncertainty of the geometric scaling and eikonalization game. 

2.4 DIS revisited 

In DIS, the photon virtuality Q 2 introduces a further scale to the process, i.e. one goes 
from Fig. |3]a to Fig. [5]a. The traditional DIS region is the strongly ordered one, Q 2 3> k\ 3> 
p 2 ± , where DGLAP-style evolution |fj] is responsible for the event structure. As above, 
ideology wants strong ordering, while real life normally is based on ordinary ordering 
Q 2 > k\ > p 2 ± . Then the parton-model description of F 2 (x,Q 2 ) in eq. (H) is a very good 
first approximation. The problems come when the ordering is no longer well-defined, i.e. 
either when the process contains several large scales or when Q 2 — > 0. In these regions, 
an F 2 (x, Q 2 ) may still be defined by eq. (0), but its physics interpretation is not obvious. 

Let us first consider a large Q 2 , where a possible classification is illustrated in Fig. ^o. 
The regions Q 2 > p\ > k\ and p\ > Q 2 > k\ correspond to non-ordered emissions, 
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that then go beyond DGLAP validity and instead have to be described by the BFKL 
17fl or CCFM || equations, see e.g. @]. Normally one expects such cross sections to be 



7 *p ^ i*emV pj T Q 2 ) - a'* p ~ rx 7 * p - rx 7 * p 4- rx 7 * p 4- rx 7 * p (22) 
°tot - /r>2 i m 2\2 2 v x )V i -°F 2 °DIS LODIS ' QCDC ' BGF " K^^J 



small at large Q 2 . The (sparsely populated) region p 2 ± > k\ > Q 2 can be viewed as the 
interactions of a resolved (anomalous) photon. 

The region k\ > Q 2 3> and k\ > p 2 ± contains the 0(a s ) corrections to the lowest- 
order (LO) DIS process 7*q — > q, namely QCD Compton 7*q — > qg and boson-gluon 
fusion 7*g — > qq. These are nothing but the direct processes 7q — > qg and 7g — ■> qq 
extended to virtual photons. The borderline k\ > Q 2 is here arbitrary — also processes 
with k\ < Q 2 could be described in this language. In the parton model, this whole class of 
events are implicitly included in F2, and are related to the logarithmic scaling violations 
of the parton distributions. The main advantage of a separation at k± = Q thus comes 
from the matching to photoproduction. Also the exclusive modeling of events, with the 
attaching of parton showers of scale Q 2 to DIS events, is then fairly natural. 

The DIS cross section thus is subdivided into 

47r 2 a em Q 2 
(Q 2 + m 2 ) 

The (T^jg is given by eq. @, while the last two terms are well-defined by an integration 
of the respective matrix element ||. When extended to small Q 2 , these two terms will 
increase in importance, and one may eventually encounter a er^Dis < 0> ^ calculated by 
a subtraction of the QCDC and BGF terms from the total DIS cross section. However, 
here we expect the correct answer not to be a negative number but an exponentially 
suppressed one, by a Sudakov form factor. This modifies the cross section: 

7*P _ 7*P 7*P 7*P 7*P I QCDC BGF 1 (0"i\ 

<T LO DIS — a DIS — ^QCDC — ^BGF ¥ a DIS eX P r^, I ' \ Z6 ) 

V ^DIS / 

Since we here are in a region where u^jg a F 2 P ^ ^- e - where the DIS cross section is 
no longer the dominant one, this change of the total DIS cross section is not essential. 
Even more, for Q 2 — > we know that the direct processes should survive whereas the 
lowest-order DIS one has to vanish. Since eq. (Dl) ensures that cr^jg — > in this limit, it 
also follows that c^odis does so. 

2.5 Prom Real to Virtual Photons 

It is now time to try to combine the different aspects of the photon, to provide an an- 
swer that smoothly interpolates between the photoproduction and DIS descriptions, in a 
physically sensible way. 

A virtual photon has a reduced probability to fluctuate into a vector meson state, and 
this state has a reduced interaction probability. This can be modeled with the traditional 
dipole factors |TB] 



-v*m P d(^ 2 )= E ^{-Jt^) a ^ 8) - (24) 
Similarly, the GVMD states are affected, 
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where a relation 2k± ~ m is assumed. 

The above generalization to virtual photons does not address the issue of longitudinal 
photons. Their interactions vanish in the limit Q 2 — > 0, but can well give a non- negligible 
contribution at finite Q 2 flnj. A common approach is to attribute the longitudinal cross 
section with an extra factor of ry = ayQ 2 /rriy relative to the transverse one ||2U|| , where 
ay is some unknown parameter to be determined from data. Such an ansatz only appears 
reasonable for moderately small Q 2 , however, so following the lines of our previous study 
of jet production by virtual photons ||, we will try the two alternatives 

2 2\ ArriyQ 2 
n(m v ,Q ) = * {< + Q2)2 , (26) 

r 2 K,Q 2 ) = a (m ? Q +Q2) • (27) 

While ri vanishes for high Q 2 , r 2 approaches the constant value a. The above VMD 
expressions are again extended to GVMD by the identification my « 2k ±. The cross 
sections can then be written as 



^vm P d(^ 2 ,Q 2 ) = £ ^ [l + r t (m 2 v ,Q 2 )} 



X 



I' 



"' T ]2 ^(W 2 ), (28) 



m 2 v + Q< 



2 If -° vm) ^ ' < 29 > 

Note that we also here replaced s by the conventional DIS variable W 2 ; obviously both 
refer to the same 7*p squared invariant mass. 

The extrapolation to Q 2 > is trivial for the direct processes, which coincide with 
the DIS QCDC and BGF processes. The matrix elements contain all the required Q 2 
dependence, with a smooth behaviour in the Q 2 — > limit. They are to be applied to the 
region k± > max(fci, Q) (and k± > p±, as usual). 

Remains the LO DIS process. It is here that one could encounter an overlap and 
thereby double-counting with the VMD and GVMD processes. Comparing Fig. |5|b with 
Fig. f|, one may note that the region p± > k± involves no problems, since we have 
made no attempt at a non-DGLAP DIS description but cover this region entirely by the 
VMD/GVMD descriptions. Also, if Q > ki, then the region k\ < k±_ < Q (and k± > p±) is 
covered by the DIS process only. So it is in the corner k± < k\ that the overlap can occur. 
If Q 2 is very small, the exponential factor in eq. ( p3l) makes the DIS contribution too small 
to worry about. Correspondingly, if Q 2 is very big, the VMD/GVMD contributions are 
too small to worry about. Furthermore, a large Q 2 implies a Sudakov factor suppression 
of a small kj_ in the DIS description. If W 2 is large, the multiple-interaction discussions 
above are relevant for the VMD/GVMD states: the likelihood of an interaction at large 
Pj_ will preempt the population of the \ow-p± region. 

In summary, it is only in the region of intermediate Q 2 and rather small W 2 that we 
have reason to worry about a significant double-counting. As we shall see, it is indeed in 
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this region that the VMD/GVMD contribution does not appear to dampen quite as fast 
as data indicates. Typically, this is the region where x ~ Q 2 / (Q 2 + W 2 ) is not close to 
zero, and where F 2 is dominated by the valence-quark contribution. The latter behaves 
roughly oc (1 — x) n , with an n of the order of 3 or 4. Therefore we will introduce a 
corresponding damping factor to the VMD/GVMD terms. The real damping might be 
somewhat different but, since small W values are not our prime interest, we rest content 
with this approximate form. 

In total, we have now arrived at our ansatz for all Q 2 

7 *p 7 *p ( °"QCDC + °"BGF \ 7 * p 7 * p 

°"tot — °blS eX P 7*p + °~QCDC + °BGF 

V ^DIS / 

+ (q^wi) T (*d + 4v P md) , (30) 

where the DIS, VMD and GVMD terms are given by eqs. (Q), and (|29|), respectively, 
and the QCDC and BGF terms by direct integration of the respective matrix elements 
for the region k± > max(ki,Q). To keep the terminology reasonably compact, also for 
the 7*7* case below, we will use res as shorthand for the resolved VMD plus GVMD 



contributions and dir as shorthand for the QCDC and BGF processes. Then eq. (30) 
simplifies to 

/ o- 7 * p \ * / W 2 \ n 
air = *l£ exp [—$=)+ a d \ p + -f— <r£> . (31) 



(Tii) \Q 2 + W 2 



2.6 Photon— photon collisions 

A generalization of the 7p picture to 77 events is obtained by noting that each of the 
two incoming photons is described by a wave function of the type given in eq. (TO). In 
total, there are therefore three times three event classes. By symmetry, the 'off-diagonal' 
combinations appear pairwise, so for real photons the number of distinct classes is only 
six. These are, cf. Fig. |6|: 

1. VMD x VMD: both photons turn into hadrons, and the processes are therefore the 
same as allowed in hadron-hadron collisions. 

2. VMD x direct: a bare photon interacts with the partons of the VMD photon. 

3. VMD x anomalous: the anomalous photon perturbatively branches into a qq pair, 
and one of these (or a daughter parton thereof) interacts with a parton from the 
VMD photon. 

4. DirectX direct: the two photons directly give a quark pair, 77 — > qq. Also lepton 
pair production is allowed, 77 — > , but will not be considered by us. 

5. DirectX anomalous: the anomalous photon perturbatively branches into a qq pair, 
and one of these (or a daughter parton thereof) directly interacts with the other 
photon. 

6. Anomalous x anomalous: both photons perturbatively branch into qq pairs, and 
subsequently one parton from each photon undergoes a hard interaction. 

Most of the above classes above are pretty much the same as allowed in 7p events, since 
the interactions of a VMD or anomalous photon and those of a proton are about the same. 
Only the direct x direct class offer a new hard subprocess. 

The main parton-level processes that occur in the above classes are: 
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d) e) f) 

Figure 6: Contributions to hard real 77 interactions: a) VMDxVMD, 
b) VMDx direct, c) VMDx anomalous, d) direct x direct, e) direct x anomalous, and 
f) anomalous x anomalous. Notation as in Fig. |2|. 

• The 'direct' processes 77 — > qq only occur in class 4. 

• The 'single-resolved' processes 7q — ► qg and 7g — > qq occur in classes 2 and 5. 

• The 'double-resolved' processes qq' — > qq' (where q' may also represent an anti- 
quark), qq — > q'q', qq — > gg, qg — > qg, gg — > qq and gg — > gg occur in classes 1, 3 
and 6. 

The classification of a generic Feynman diagram by the different possible components is 
illustrated in Fig. [7|. The appearance of more scales makes it infeasible to draw diagrams 
like Fig. |^b, but the principles are the same. 

Also the extension to virtual photons follows from the 7*p formalism above, but now 
with (up to) five scales to keep track of: p±, k±±, k± 2 , Qi and Q 2 . First consider the 
three by three classes present already for real photons, which remain nine distinct ones 
for Q\ 7^ Q\. Each VMD or GVMD state is associated with its dipole damping factor 
and its correction factor for the longitudinal contribution. The QCDC and BGF matrix 
elements involving one direct photon on a VMD or a GVMD state explicitly contain the 
dependence on the direct photon virtuality, separately given for the transverse and the 
longitudinal contributions. Also the direct x direct matrix elements are known for the four 
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7 W\A/ 



7WVV 



^±2 



1. VMDxVMD: k ±i ,k ±2 < h, arbitrary p ± 

2. VMDx direct: k±i < k < p± < k ±2 
+ (1 <->2) 

3. VMD x anomalous: k±i < ko < k± 2 < P± 
+ (1 <->2) 

4. DirectX direct: k < k±i = k± 2 

5. DirectX anomalous: k Q < k^i < p± < k±_ 2 
+ (1 <->2) 

6. Anomalous x anomalous: k Q < k±i,k± 2 < p± 



Figure 7: Schematic graph for a hard 77 process, showing the three different scales. 
To the right is shown the relation to the six classes in the text. 



possible transverse/longitudinal combinations. Some examples should be enough: 



^Jdxgvmd^ 2 , Ql Ql) = E [1 + g?)] 



2 \ 2 

'v, 



^IX^f [l + r,(4*l,«S 



x i«rfoi) 2 % 5I ' T> '' l ' lM>(w ' 2) ' (32) 



2 \ 2 



X 



x {a^{W\Ql) + al^{W\QD) . (33) 

To this should be added the new DIS processes that appear for non- vanishing Q 2 , 
when one photon is direct and the other resolved, i.e. VMD or GVMD. For simplicity, first 
assume that one of the two photons is real, Q\ = 0. For large Ql, this DIS contribution 
can be given a parton-model interpretation, 

(34) 

Note that this is only the resolved part of (T^ot ■ The direct contribution from 7*7 — ► qq 
comes in addition, but can be neglected in the leading-order definition of F%. We will 
therefore use parton distribution parameterizations for the resolved photon, like SaS ID 
Tjfl , to define the c r DiSxres( ( 5i)- Then eq. (|3l|) generalizes to 



a DISxres , 



/ W 2 \ n 

~ u resxdir ~ "dirxdir ~ I q 2 _|_ yy 2 J "resxres ■ \' J ' J J 
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The large-x behaviour of a resolved photon does not agree with that of the proton, but 
for simplicity we will stay with n — 3. 

The generalization to both photons virtual then gives 



where the choice of damping factor for the last term again is a simple guess for an exten- 
sion. When Ql S> Q 2 the expression for cr^ 7 (W 2 , Q\, Q 2 ,) can be related to the structure 
function of a virtual photon, F 2 7 *(x, Q 2 = Q 2 , P 2 = Ql), where x = Q\/{Q\ + Q\ + W 2 ). 

A comment about the exponential factors suppressing the DIS terms. Properly, each 
VMD/GVMD state should come with its ratio of direct to DIS cross sections. However, a 
number of common factors divide out in this ratio: obviously the probability to fluctuate 
to the state in question, but partly also the form of the parton distributions in the 
state. Therefore, it is a good approximation to define a common exponential form for 
all VMD/GVMD state, based on the weighted average, which is equivalent to using the 
full resolved term in both numerator and denominator. 



7 7 

^DISxres eX P 
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res X res ) 



(36) 



2.7 Other Model Aspects 

In the preceding section we have mainly described the model as differential in W (invariant 
mass of the 7*p or 7*7* system) and the hardest scales involved: p±, k ±i and Q 2 . The 
discussion about the different processes has been limited to the parton-level only, but it is 
necessary to add parton showers, beam remnants, hadronization etc. in order to simulate 



complete events. The model is implemented in the Pythia 6.151 event generator |21 
which includes the above and many other aspects, and it will be used in the following 
for the results presented. Clearly, here we can only give a few typical examples; with 
the help of the generator it is possible to study any specific experimental conditions and 
observables. 

While based on the formalism outlined above, the complete event simulation also 
involves some more model-specific assumptions. We here present a few of those, to try 
to complement the overall picture. Further refinements are possible in many places, and 
some of the known shortcomings are mentioned. 

The extension from pp to 7*p and 7*7* collisions requires two new aspects to be 
introduced in the event simulation structure. One is the need to construct the e — > ej* 
kinematics, and to let the collision description depend not only on the resulting W but 
also on the selected Q 2 virtuality scale(s) of the photon(s). This part was described in ||. 

The other is the necessity to mix all the different reaction processes. The main clas- 
sification is here into 4 components for 7*p — VMD, GVMD, direct and DIS — and 13 
for 7*7* — 4-squared except that there is no DIS process on a direct photon or another 
DIS photon. Each of these components has a set of allowed subprocesses, e.g. 7*q — > qg 
and 7*g — > qq for direct processes. This set is the same for a VMD and a GVMD photon, 
but the two are distinguished by the different parton distributions and total cross section 
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parameterizations used. In the Monte Carlo, a choice is first made between the 4/13 
components and then, for that component, among the allowed subprocesses, according to 
maxima estimated at the initialization. Once the kinematics of the event has been fully 
selected, the ratio of the actual cross section to the assumed one is used to retain the 
event or to select a new component and subprocess. 

The separation of VMD from GVMD requires access to parton distributions where 
those two components are explicitly made available separately. Effectively this limits the 



choice to the SaS parameterizations and with the assumed k « 0.5 GeV to SaS 
ID. An almost equivalent formalism could have been constructed in terms of a common 
resolved class, so that we could have gone from 4/13 to 3/6 categories. While the cal- 
culation of partonic processes would have been simplified, aspects related to total cross 
section processes (elastic, diffractive, \ow-p±, etc.) would not, so we have not yet tried to 
construct such a complete alternative machinery. 

A VMD photon may be classified either as p°, u, <f> or J /if). The parameterizations 
of total, elastic, single and double diffractive cross sections for V^p and V1V2 collisions 



of real photons are given in [fL5|| . The respective inelastic non-diffractive cross section 
is obtained by subtraction, and sets the envelope within which the jet cross section is 
eikonalized. The events are ultimately classified, either by the hardest interaction that 
occurs, or as a low-p_i_ event if there are no hard interactions. For virtual photons, both 
the total, elastic, single and double diffractive cross sections are dampened by the same 
dipole factors, cf. eq. (|24j). The jet cross section is obtained by the virtuality-dependent 
parton distributions, however, so in the eikonalization procedure the mixture between jet 
and low-p_i_ events is explicitly Q 2 -dependent. For technical reasons, currently it is not 
possible to mix different dipole factors for the different VMD states, so the p/uo mass is 
used throughout, thus giving a too fast dampening of the (rare) J /if) events. 

For the GVMD states, the mass selection is based on the identification m ~ 2k± and a 
k± spectrum e.g. as implied by eq. (|29|) . Thus, neglecting the longitudinal-photon factor, 
one has a mass spectrum like dm 2 /(m 2 + Q 2 ) 2 in the range 2fc < m < 2p^ min (W 2 ). This 
spectrum begins at 1 GeV, i.e. a bit above m p . Again there is not yet any provision for 
heavier quarks in the dipole dampening formulae, but the mass of an ss state is shifted 
by rn^ — m p and that of an cc by mjw — m p , to ensure that the GVMD spectrum starts 
above the lowest-lying state. The mass selected above then becomes the final-state mass 
of an elastically scattered GVMD particle. A diffractively scattered GVMD is assigned 
an excited mass according to a spectrum stretching from the selected mass and upwards 
to the kinematical limit, essentially like dm 2 /m 2 , as described in fTB] . 



An elastic or diffractive GVMD state is in reality no different from a diffractive VMD 
state of the same mass. A low-mass system is allowed to decay to two hadrons, whereas 
more massive ones are considered as strings stretched along the event axis, modulo limited 
transverse-momentum fluctuations. It is assumed that a simple string is stretched between 
two endpoint quarks half of the time, and a hairpin string arrangement with a gluon pulled 
back the other half. Clearly this mixture is only a simple first approximation to what is 
likely to be a more complex structure [2l| . 

The mass of the GVMD state does not enter the description of inelastic non-diffractive 
events. However, here the k±_ ~ m/2 does provide a 'primordial k± kick that can be 
transmitted to the parton of a hard scattering, not only for GVMDxp but also in processes 
such as DISxGVMD. The corresponding number ought to be k±_ m m p /2 for VMD and 
k± « m p /3 ~ m p /2 for a proton. However, in hadronic collisions much higher numbers 



than that are often required to describe data, typically of the order of 1 GeV |23[ if a 
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Gaussian parameterization is used. Thus, an interpretation as a purely nonperturbative 
motion inside a hadron is difficult to maintain. 

Instead a likely culprit is the initial-state shower algorithm. This is set up to cover 
the region of hard emissions, but may miss out on some of the softer activity, which 
inherently borders on nonperturbative physics. By default, the shower does not evolve 
down to scales below Qo — 1 GeV; if the kj_ ~ m/2 scale of a GVMD photon is above 
this the evolution is stopped already at the larger scale. (Whether Q Q should best be 
compared with m or m/2 is an open issue.) Any shortfall in shower activity around or 
below this cutoff then has to be compensated by the primordial k± source, which thereby 
largely loses its original meaning. 

It could be argued that, while a VMD photon should have the same primordial kj_ 
spectrum as a hadron, a GVMD one should receive its main contribution from the per- 
turbative k± ~ m/2 contribution, i.e. have a shape like dk\/(k"j_ + Q 2 /4) 2 in the range 
k < k± < p± m m(W 2 ). However, this has the questionable consequence that a low-mass 
GVMD state would have a smaller average 'primordial k±' than a VMD one. While we 
have retained such a possibility as an option, our default has instead been to go to the 
other extreme, where the partons in a GVMD state has a k± given by the vector sum 
of the perturbative power spectrum and the same nonperturbative Gaussian smearing as 
for hadrons. Comparisons with data on photon remnant jets should eventually shed more 
light on the appropriate procedure. 

In most processes, initial- and final-state shower activity is routinely added to the 
lowest-order process, thereby providing an approximation to higher-order QCD correc- 
tions. An exception is the DIS process 7*q — > q, where currently only final-state radiation 
has been implemented. The technical reason why the initial-state radiation algorithm does 
not work here is that it is based on a definition of the z of the splitting kernel as being 
the fraction by which the s of the hard subprocess is reduced by an emission, and s = 
in the DIS process above. The traditional z interpretation is sensible if one remembers 
that the complete DIS hard process is eq — > eq rather than only 7*q — > q. However, since 
we have split off the photon flux, the complete process is not trivially available currently. 
For small Q 2 values, this is no major problem, since events with activity above this scale 
belong to the direct event class. Therefore the shower is constrained only to populate the 
region below Q 2 . The problem may become more severe in the high-Q 2 domain, where 
DIS processes provide the dominant event class and there is a large phase space for shower 
activity. That is not the region of main interest to us here, however. 

The choice of a formalism, with the virtual photon flux separated from the hard pro- 
cesses, is there for two reasons. One is that it makes it possible to set up selection criteria 
for an event sample, e.g. on W and Q 2 or in terms of scattered electron energies and 
angles, that are consistent across the different photon components. Another is that many 
nonperturbative physics aspects, such as VMD total cross sections, best are formulated 
in such terms. 

Also a few other areas are not fully developed. One example is multiple interactions, 
where only the simpler impact-parameter- independent option can be used J7|, at least 
until one has constructed a model of the spatial distribution of partons in a GVMD state. 
On the positive side, options are available for several of the other aspects discussed in 
this article, so that one may study the sensitivity to the assumptions made. 
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Figure 8: a) The total <r7 ot p cross section as function of the photon virtuality at 
Wry* p = 100 GeV. The resolved components shown are with alternative r\ for the 
longitudinal contribution, b) Contributions to the total o"7 ot p cross section from the 
different event classes: direct, VMD, GVMD and DIS. For simplicity, the VMD and 
VMD+GVMD contributions are not shown for the lines corresponding to Q 2 values 
in the range 0.85 — 2.5 GeV 2 . See text for details. The legend for the data points 
can be found in Fig. |]a. 

3 Results 

In this section, some examples of event properties will be given, with emphasis on the 
differences between the direct, VMD and anomalous/GVMD component of the photon. 
We will concentrate on a few distributions, starting in the photoproduction limit, and then 
explore how they gradually turn over from being dominated by the resolved processes to 
finally being dominated by the bare direct and DIS ones as the photon virtuality increases. 

Owing to the resolved components of the photon, pp, 7p and 77 events show similar 
behaviour for several event properties, such as multiplicity and transverse energy flow. 
But clearly, differences are expected due to, for example, the additional direct event class 
in 77. 

3.1 Total Cross Sections 
3.1.1 7*p Total Cross Section 

The SaS ID photon parton distribution and the CTEQ5L proton parton distribution will 
be used throughout in this section if not otherwise stated JTT|, p4[ . In Fig. ||a, the 7*p 
cross section at W 7 * p = 100 GeV is shown as a function of the photon virtuality Q 2 . The 
VMD and GVMD components contain contributions from resolved longitudinal photons 
with the ri alternative. An implicit y dependence enters the factors, eq. (p6|) and fl27D , 
to compensate for the difference in photon flux between longitudinally and transversely 
polarized photons. To estimate the longitudinal contribution to the resolved component 
in 7*p cross sections, a fixed y is chosen in order to set the relative amount of transverse 
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and longitudinal photons. The value of y consequently depends on the beam energy, and 
at HERA W 7 * p = 100 GeV would correspond to y ~ 0.1, a value which will be used in 
the following. The parameter a, appearing in r i; is unknown and need to be determined 
by data. As a starting point for the discussion, a is chosen to 0.5. 

In the photoproduction limit, the VMD component dominates the total cross section 
whereas the DIS one has a vanishing cross section by construction, see Fig. ^|a. The 
GVMD and direct processes share the remaining part, about 20%, of the total cross 
section. It follows from the introduction of the dipole factors in eq. ([24]) that, when 
considering total 7*p cross sections, the VMD class is the dominating component up to 
photon virtualities of the order of a vector meson mass. Similarly, the characteristic 
virtualities for the GVMD class to still be important are of the order of the ki parameter 



discussed in section O 



The direct processes are simulated in the region where the hard scale of the 
photon-parton scattering fulfill kj_ > max(fci,Q), cf. Fig. ^. At low Q 2 values, the 
k\ = p±mm{W 2 ) ~ 1.3 GeV provides a fix limit on the available phase space, but at higher 
Q 2 the phase space shrinks, explaining the kink in Fig. ||a for the direct cross section at 
Q 2 = k\. So it is rather this condition that drives the dampening at large Q 2 than the 
explicit Q 2 dependence appearing in the matrix elements. A kink can also be seen in the 
DIS cross section; it is an artifact from the freezing of the parton distributions below their 
range of applicability (as introduced in section [2.1|) . Hence, the CTEQ5L proton parton 
distribution has Q 2 nin = 1 GeV. The factors introduced for the DIS component, eq. (|9]), 
together with the exponential dampening, eq. (|23|) , make the unusual form of the DIS 
cross section — to be discussed in more detail below. 

The proton structure function F<i has been measured at various x and Q 2 values [p5[| 
and is related to the total virtual photon-proton cross section through eq. (|^). In 
Fig. ||b and|9]a, o"7 ot p is shown from zero Q 2 , the total photoproduction cross section p6| , to 



medium Q 2 values as a function of the invariant mass W of the 7*p-system. Results from 
the model are compared with data from different fixed-target and HERA experiments. 
The contribution from the different event classes: direct, VMD, GVMD and DIS are 
shown in Fig. |]b, and sum up to the total contribution as given in eq. (|31|). The DIS 
contribution can be obtained by subtracting the 'VMD+GVMD+Dir' contribution from 
the 'Total rj.' one. With increasing photon virtuality, the energy dependence of k\ makes 
the GVMD states to be less dampened at large energies as compared to low ones. This is 
also indicated in Fig. |||b by comparing the 'VMD' and 'VMD+GVMD' lines at different 
Q 2 and W. 

In Fig. |)a, the total contribution obtained with different assumptions of the longitu- 
dinal contribution for resolved photons are compared to the case with transverse photons 
only. With the other model assumptions made, a non-vanishing longitudinal contribution 
is indicated. The alternative with r x gives a good fit to data, but overshoots at low W 
values for some Q 2 values. In the high Q 2 region, showed in Fig. §a and at the bottom 
in Fig. [5|b, the choice of longitudinal alternative for the resolved component is irrelevant 
for the total cross section since the resolved event classes give a negligible contribution. 
Hence, the parameter a should be tuned at intermediate Q 2 values. At lower Q 2 values, 
however, the r 2 alternative is overshooting data significantly. Taking a = 0.5 as a satis- 
factory value for the r\ alternative, the r 2 one require a = 0.2 to give an equally good 
description (not shown). 

The direct processes increase with W 7 * p faster than the DIS one, and consequently 
also the exponential suppression of the DIS term. Therefore they contribute substantially 
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Figure 9: a) The total contribution to o"7 ot P from the different event classes. The 
two different longitudinal alternatives are compared to the case with only transverse 
photons considered. The data and 'Total r\ is the same as in Fig. ||b. b) Same as in 
a but for a different Q 2 range and also the DIS contribution is shown. For reference, 
the Q 2 = 3.5 GeV 2 lines are the same, however. 

to the total cross section in the large W 7 * p (small x) region, even at photon virtualities 
up to tens of GeV 2 , shown for W 7 * p = 100 GeV in Fig. ||a. The DIS and direct classes 
are often combined and clearly dominate the total cross section at Q 2 = 3.5 GeV 2 which 
is above the my and k\ scales of the VMD and GVMD states discussed above. Due to 
the negligible contribution from the direct processes at W 7 * p = 10 GeV, the exponential 
suppression factor for the DIS term is close to unity — independently of Q 2 . Hence, the 
sum of the direct and DIS contributions reproduce the a^jg cross section, eq. (]22|) . For 
W r * p = 100 and 1000 GeV this is the case for Q 2 > 5 and 8 GeV 2 , respectively. 

3.1.2 7*7 Total Cross Section 

For simplicity, a fixed yi — 0.1 is also used in photon-photon collisions for the calculation 
of the longitudinal resolved photon contributions. In Fig. pH|a, the a2 Q t at W 7 * 7 = 100 GeV 
is shown as a function of one of the photon virtualities; the other photon is real. The 
different components are shown separately and sum up to the total contribution with 
the ri alternative. Additionally, relative to the corresponding 7*p case, the direct cross 
section enters and starts to be significantly dampened only when Q\ ~ W 2 . The kink 
seen for the single-resolved case can again be explained by the constraint of the hard scale 
in the scattering process to be larger than the virtuality of the direct photon. 

The same distribution is shown for the case of two virtual photons in Fig. |10|b>. The 
low virtuality end of Q\ of course corresponds to the cross sections at Q\ = 10 GeV 2 in 
Fig. |l0la. The DIS cross section dominates in this region, with the first photon as the 
target probed by the second one. This is interchanged for Q\ > Q\. It is noticeable that 
all the different processes shown are of about the same importance when Q\ is in the 
neighbourhood of Q\. 
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Figure 10: a) The cr^ 7 at W 1 * 1 = 100 GeV as a function of the photon virtuality Q\ 
of one of the photons is shown for the DIS, direct, single- resolved (1-res), double- 
resolved (2-res) and total contribution. The other photon is real, Q\ = GeV 2 . b) 
Same as in a) but for 7*7* with a target photon virtuality Q\ = 10 GeV 2 . 



For various reasons, the photon K7 measurements |£7 are both less precise and avail- 
able in a smaller kinematical range than in the proton case. In Fig. |TT|a, the is shown 
as a function of W for different Q 2 values, similarly to the study in the previous section. 
With a = 0.5 for both alternatives of the longitudinally resolved photons, data does not 
discriminate between them. The r\ alternative is below the high-energy end data points 
for the Q 2 = 4.04 and 5.15 GeV 2 lines but the errors and spread in data do not give an 
unambiguous conclusion. 

for some different event classes at Q 2 = 



lib shows the contribution to a 



1 7 



Fig 

9.90 GeV 2 '. The direct and DIS ones are dominating at low 7*7 invariant masses whereas 
the single-resolved ones increase in importance at high W 1 * 1 . The direct processes show 
the characteristic 1/W 2 fall off || and thereby give its major contribution to the region 
where the valence quarks dominate the contribution to F^. Again, the single-resolved 
processes (1-res; cf. direct ones in 7*p) increase with energy, which therefore dampens 
the DIS term slightly at large W 7 * 7 , a region which is dominated by the gluon content 
of the photon. As resolved processes in 7*p, double-resolved ones (2-res) in 7*7 increase 
with Wry* 7 . 



3.1.3 7*7* Total Cross Section 

There are hardly any systematic data on 7*7* cross sections. Prior to the recent and ongo- 



ing LEP studies, essentially the only publication is the F e g measurement by PLUTO |28| 



A comparison with these data points is found in Fig. |12|a. The low W values — the first 
point is in the resonance region — makes the comparison especially precarious, and the 
second photon is also not all that virtual with its = 0.35 GeV 2 . The main point 
of the plot is thus not the acceptable agreement with data, but to illustrate the amount 
of reduction of the cross section relative to the real-photon alternative Q\ = and the 
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Figure 11: a) Same as in Fig. |^b but for 7*7 instead of 7*p. b) The a^ t as a function 
of the invariant mass of the 7*7 system is shown for the DIS, direct, single-resolved 
(1-res), double-resolved (2-res) and total contribution. 



dependence on assumed longitudinal contributions. 

Double-tag two-photon events, e + e~ — > e + e~ + hadrons, have been measured by the 
L3 collaboration p9l- The differential e + e~ cross section w.r.t. the variable 

y = i« I I (37) 



■MS 



is shown in Fig. Ob. The two scattered leptons are required to be within the polar angle 



of 30 < 9i < 60 mrad (with respect to the incoming beams) and to have an energy larger 
than 30 GeV. This plot is expected to provide a clean test of the BFKL behaviour, with a 
cross section increasing at large Y ||30| , but the data does not support such a behaviour, 



or at least less of it than expected. Our model agrees with the data at small Y, but 
tends to fall below at larger values. This is almost entirely dictated by the drop in the 
dominant direct contribution. The resolved contributions — which are the ones that could 
be used to represent a BFKL behaviour in our framework — do come up at large Y, but 
nowhere near enough to make a significant contribution. We remind that the resolved 
contributions are suppressed in Q 2 by simple dipole factors; a less steep drop and thereby 
a larger resolved contribution could be motivated e.g. if the eikonalization of the jet cross 
section is less significant for a virtual photon. Instead of the conventional increase of the 
VMD cross section like s e with e ~ 0.08, eq. (0), a larger effective value like e ~ 0.2 — 0.3 
could thus easily be accommodated for virtual photons pTJ. For the moment, however, we 
prefer to await further experimental data. In particular, more detailed analyses of event 
properties could provide some insight on the relative mixture of event classes at large Y. 

3.2 Event Properties 

Some event properties are fairly similar between pp, 7*p and 7*7*, while others differ 



markedly. The charged multiplicity, Fig. 13, is an example of the former, while the jet 
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Figure 12: a) The oy ot 7 cross section as a function of W^*^* for a probing and target 
photon with virtualities 5 and 0.35 GeV 2 respectively. The data points from PLUTO 
are obtained from the measured F e g through eq. (|7|). For reference, the results with 
a real photon target is shown, Q\ = b) The differential e + e~ — > e + e~ + hadrons 
cross section as a function of Y — \n{y\y2S / Q\Q%) . The different components 
shown add up to the 'Total n'. 



rate, Fig. [14], is of the latter. 

Over a wide range of Q 2 values, the 7*p and 7*7* events have average charged multi- 
plicities within ±20% of the pp value, Fig. |13|a. (Elastic and diffractive events have not 
been considered. Their inclusion would reduce average multiplicities and change some 
details, but leave a similar overall picture.) For 7*p at small Q 2 , this is a consequence 
of the dominance of VMD events, with the characteristic features of hadronic physics, 
including such aspects as multiple parton-parton interactions. The VMD multiplicity is 
fairly Q 2 -independent, Fig. pT3|b. The drop to a somewhat lower 7*p multiplicity around 
Q 2 pa 1 — 10 GeV 2 instead comes from the transition from VMD dominance to DIS ditto, 
where the DIS events are cleaner by consisting of only one string piece stretched directly 
between the kicked-out quark and the beam remnant. Remember that the DIS events 
here by definition are low-p± ones, with the high-p^ part in the direct class, with a higher 
multiplicity. As the borderline between the two shifts with Q 2 , both individually show 
an increasing multiplicity: the direct by corresponding to a smaller fraction of higher-px 
events and the DIS by allowing an increasing admixture of jet events. Figures |T3|c and |l3|d 
show the corresponding Ql dependence of the 7*7* multiplicity, for a fixed Q| = 1 GeV 2 , 
normalized to the number of event of each kind and to the total number, respectively. 
Again the disappearance of the doubly-resolved event class is responsible for a drop in 
the multiplicity, with the DIS and direct processes taking over. We remind that in 7*7* 
the high-px part of the DIS process is the single-resolved one, corresponding to the direct 
one of 7*p, while the 7*7* direct class has no correspondence in 7*p. 

When it comes to the fraction of events that contain at least one jet, the differences 
are orders of magnitude, both with a hierarchy of increasing jet fraction from pp to 7*p 
to 7*7*, and with an increase as a function of Q 2 , Fig. |l4|a. Again the detailed studies in 
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Figure 13: a) The average charged multiplicity (n c h) as a function of the photon 
virtuality for 7*p, 7*7 and 7*7* events at a center of mass energy of 100 GeV. The 
result from pp events is indicated with a cross on the y-axis. For the 7*7* events, 
the other photon virtuality is kept fixed at Q\ = 1 GeV 2 . (In 7*7 Ql = GeV 2 .) 
No diffractive or elastic events are considered, b) The results from the different 
components in 7*p, averaged over the number of events of the respective kind, c) 
The result from different components in 7*7*, averaged over the number of events of 
the respective kind, d) As in c) but averaged over the total number of 7*7* events. 
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Figure 14: a) The fraction of events that contain at least one jet with an E± > 5 GeV 
inside a cone of radius R = sj (Ay) 2 + (A0) 2 < 1. W = 100 GeV in all cases. 
No diffractive or elastic events are considered, b) The results from the different 
components in 7*p, averaged over the number of events of the respective kind, c) 
The result from different components in 7*7*, averaged over the number of events of 
the respective kind, d) As in c) but averaged over the total number of 7*7* events. 
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Figs. |H|b-d show that the main effect in Q 2 is not in each individual event class but rather 
in the mixture of the event classes. Thus VMD is similar to pp, with minor differences 
from parton distributions, while GVMD is higher and direct higher still (at large Q 2 ). The 
DIS rate starts vanishing, since events with parton p± > Q are put in the direct class, 
but eventually comes up for Q > p±. In that range, all direct events are required to have 
jets at lowest order, and only fluctuations in the parton shower and hadronization can 
occasionally stop jets from being found. For 7*7* events, shown in Figs. 14c and d, the 
single- and double-resolved processes show the behaviour of the direct and resolved (VMD 
plus GVMD) ones, respectively, in 7*p. The additional direct class increase slowly with 
Q 2 , Fig. |T4]c, and dominates the jet rate between the single-resolved and DIS dominated 
regions, peaking at Q\ = 1000 GeV 2 , Fig. |T4]d. It is at this point, for this center of mass 
energy, where the direct matrix elements starts to be suppressed by the large photon 
virtuality. 

The transverse energy flow, as a function of pseudo-rapidity in the center of mass 
frame of the collision, is shown in Fig. [15|. (Again, elastic and diffractive events are not 
considered.) In Fig. [l5| a all photons are real except for 77* where one of them has the 
virtuality Q\ = \ GeV 2 . Here the most interesting aspect is to see how the events differ 
in the photon and the proton directions, pp being symmetric as well as 77 whereas the 
7P interpolate between the two with the photon in the rj > direction. Comparing 
the 77* with the 77 case, the single-resolved processes in 77* has a larger energy flow 
in the 7* direction (rj < 0) due to the increased jet activity. On the other hand, for 
the double-resolved ones it decreases in both hemispheres but less in the 7* direction. 
Additionally, the DIS processes come in to play but add on more or less evenly between 
the two hemispheres at this low virtuality. 

In Fig. |H)|b, the transverse energy flow is shown for three different virtualities in 
the 7*p center of mass frame. The major differences again come from the transition 
of the events from being VMD dominated to being dominated by the DIS processes. 
This is compensated by an increasing jet activity in the photon direction, giving the 
rather subtle changes of the total contribution in the photon hemisphere which camouflage 
the big changes of the relative composition of the different event classes. The different 
components in 7*p collisions with Q\ = 4 GeV 2 are shown in Fig. |l5|c, normalized to the 
number of events of the respective kind. The asymmetry is largest in the direct class, 
where all of the photon energy goes into the two high-p^ jets. The typical p± of the 
jets is not all that high, however, so at least one of the jets appears at large positive 
pseudo-rapidities. The VMD and GVMD events are more symmetric, but the virtual 
photon parton distributions are harder than the corresponding proton one, giving rise 
to the observed asymmetry. The VMD events are nearly symmetric for the real photon 
case, whereas for the GVMD ones the asymmetry is still pronounced, as it should be. As 
expected, the largest contribution in the photon hemisphere is in decreasing order from 
DIS, direct, GVMD and VMD events, where the direct and GVMD events are of about 
equal importance. 

A slightly asymmetric energy sharing between the two coloured beam remnants in a 
hadron-hadron kind of collision has been chosen for the results in this section. A more 
asymmetric energy sharing could also well be imagined, leaving an uncertainty in the 
model. For example, the energy flow in the VMD events of Fig. |l5| c could well decrease 
by several percent for a latter case. 

In Fig. |l5|d, the transverse energy flow for the collision of two photons is shown. One 
of the photons is kept at Q 2 , = 1 GeV 2 (77 < 0) and the other has Q\ = 0, 1 and 4 GeV 2 , 
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Figure 15: a) The transverse energy flow as a function of pseudo-rapidity in the 
center of mass frame of the collision, W = 100 GeV in all cases. No diffractive 
or elastic events are considered, b) The transverse energy flow as a function of 
pseudo-rapidity in the 7*p center of mass frame, normalized to the number of events. 
Distributions are shown for photon virtualities of 0, 1.0 and 4.0 GeV 2 . c) The 
different components in 7*p collisions, normalized to the number of events of the 
respective kind. The photon virtuality is Q\ = 4 GeV 2 . d) The transverse energy 
flow as a function of pseudo-rapidity in the 7*7* center of mass frame, normalized 
to the number of events. Distributions are shown for photon virtualities Q\ of 0, 1.0 
and 4.0 GeV 2 , the other photon has the virtuality Q\ = 1 GeV 2 . 
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Figure 16: a) The transverse momentum spectra of charged particles for different 
collisions, normalized to the total number of events. The pseudo-rapidity interval is 
constrained to [0.5,1.5] (corresponding to the 7 direction in case of 7p, 77 and 77* 
collisions). For 77*, one photon is real and the other has a virtuality of Q\ = 1 GeV 2 . 
b) With Q\ = 4 GeV 2 and Q\ = \ GeV 2 , the 7*7* distribution is shown for the 
different components; single-resolved, double- resolved, DIS and direct. The sum of 
them, add up to the full line. 



all shown separately. The Q\ = GeV 2 distribution was discussed earlier in compari- 
son with two real photons, and is built up by (in order of importance) the asymmetric 
double- and single-resolved contributions together with a small symmetric DIS contribu- 
tion — the direct event class is negligible. With both photons having the same virtuality, 
Q\ = 1 GeV 2 , the above relative importance between the different components still holds 
but now all event classes give a symmetric contribution to the E± flow. The last case, 
with two photons at different virtualities, the single-resolved processes are responsible for 
most of the asymmetric shape of the total contribution. Now also the direct processes 
starts to be important with a central plateau of two units in pseudo-rapidity. It is here 
comparable with the DIS contribution, about 0.25 GeV at central rapidities (normalized 
to the total number of events). 

The transverse momentum spectra of charged particles become harder when going 
from pp to 7p to 77*, Fig. [H| which is partly a reflection of the property of the respective 
parton distributions and the pattern seen in the jet rates. In 7p, the VMD processes 
dominate in the low end of the spectra and the GVMD and direct processes in the high-p^ 
tail. With increasing photon virtuality, the DIS processes enters but is only of importance 
at low p± due to the constraint p± < Q. (The parton shower and hadronization will cause 
some particles to be found at p± > Q, however.) With a photon virtuality of a few 
GeV 2 , a similar spectra is obtained as for 7*7* in Fig. |16)d, remembering to associate the 
7P direct events with the 7*7* single-resolved ones, and the 7p resolved with the 7*7* 
double-resolved ones. 

With two photons, the direct processes give a tail out at high p± that is comparable 
to the single- and double-resolved processes, Fig. |T6|p, and contribute to the hardening 
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of the spectra as compared to pp and 7p collisions. When both photons have the same 
virtuality, for example two real photons, the DIS processes are as usual absent. Although 
the mixture of event classes differs significantly within a certain kind of collision, the low 
end of the spectra remain approximately the same when the photon virtuality is increased. 



4 Summary and Outlook 

We have in this article tried to outline a scenario that covers 'all' photon interactions, 
whether real or spacelike virtual, that produces hadronic final states. Part of it relies on 
previous studies, e.g. on jet production with virtual photons or total cross sections for real 
photons, but here it is eventually combined to give the overall picture. In doing this, we 
attempt to integrate various aspects into a reasonably consistent overall picture. It would 
have been nice if a single very economical ansatz for the photon could be made to cover 
all relevant phenomenology. We do not exclude that this would be possible, although we 
could also see problems with such an approach. For instance, it appears likely that there 
really is a fundamental separation into a 'resolved' and a 'direct' (or 'pointlike') part of the 
photon, e.g. based on the x 7 distributions observed at HERA []32| |. Even this separation, 



of course, is only theoretically well-defined to lowest order, and needs to be prescribed in 
higher orders of perturbation theory. Nevertheless, the 'all orders' data quite nicely show 
a separation. 

Our model takes such an approach one step further. The resolved part of the photon is 
subdivided into VMD and GVMD, depending on whether the resolved photon is associated 
with one of the lowest-lying vector mesons or with one in the set of not so well known 
higher-mass resonances. For a virtual photon, these states then have a dampened cross 
section given by dipole form factors. The direct sector is somewhat more complicated, 
since the lowest-order DIS process 7*q — > q is not allowed in the limit Q 2 —>■ while the 
higher-order direct ones 7*q — > qg and 7*g — > qq are. It therefore requires some care to 
retain only the latter for Q 2 = while equating them with the first-order QCD corrections 
to the DIS process for large Q 2 . 

The first test that the mixing makes sense comes from a comparison with the total 
cross sections of 7*p and 7*7* interactions. It is possible to obtain a reasonable descrip- 
tion of all the data, although with some disagreement for rather small W, where our 
language is not expected to survive limited phase-space corrections and exclusive final- 
state effects anyway. Of course, since we have not even attempted to produce our own 
sets of parton distributions, but taken existing ones, all the possibilities of tuning are 
not exhausted by far. Also the inclusion of the effects of longitudinal resolved photons 
remains an area where little is known, and our simpleminded ansatz could be made more 
sophisticated, e.g. by the use of parton distributions designed to describe the partonic 
structure of longitudinal virtual photons p3| . However, the model is so far limited to the 



SaS ID set of parton distributions. An alternative approach with a common resolved class 
could have been taken, facilitating the handling of the partonic processes, but probably 
ending up with a more complicated picture for other aspects such as elastic and diffrac- 
tive scattering, multiple parton-parton interactions, etc. and at the same time losing the 
flexibility of having the non-perturbative part of the photon structure separated from the 
perturbatively calculable one. It is a subject for further studies, however. 

More sophisticated tests come from the study of event shapes. We have in this article 
provided some examples how event properties vary between different initial states and 
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photon virtualities. The examples are mainly chosen as simple illustrations; we look 
forward to more detailed studies by the experimental community, based on the code we 
now provide in the Pythia event generator. Again one cannot expect perfect agreement, 
but at least an overall such, where disagreements hopefully could help provide hints in 
which direction to move for an even better and more complete picture. In this sense, our 
model could be seen as a straw man, although a rather more sophisticated such than is 
often the case. 

Given the fairly complex description, we do not expect the model to be competitive 
in the high-Q 2 region of HERA, say for Q 2 > 10 GeV 2 . There the simple DIS language 
provides a powerful starting point, that has been well studied over the years, with many 



models developed in extensive detail |34J|. We are here more interested in the crossover 
region, Q 2 ~ 0.5 — 5 GeV 2 , where our model predicts all of the photon components to 
have comparable cross sections, and we therefore would expect no simple picture to work. 
In 7*7* studies at LEP2, it is enough to have one photon in this region for our approach 
to offer interesting alternatives to other descriptions |35| . For smaller Q 2 , it smoothly 
attaches to the existing Pythia model of real-photon interactions, while the high-Q 2 end 
has no such correspondence. 

There are some areas where we already now know that not enough effort has gone 
in to cover the field. One such is the treatment of heavy flavours and in particular 
charm production, both open and closed (primarily J/ip), where mass effects are very 
important. For the lighter quarks, like the u one, the intrinsic 'current algebra' mass is 
negligible relative to the 'constituent' one, that sets the scale e.g. of the p mass. But 
for a complete charm description, both the impact of the current algebra mass scale and 
of the further QCD-induced confinement mass effects have to be considered. Another 
example of a missing area is that related to rapidity gaps, which are already included in 
the VMD/GVMD sector within a traditional Regge language, but currently not in the DIS 
region, in spite of the quite conclusive evidence for this from HERA |36j . The two aspects 
come together in the W dependence of the production rate of exclusive vector mesons, 
where data show a steeper rise for J/i(j than lighter mesons, and a steeper rise also at 
larger Q 2 [[H]]. On the technical side, initial-state radiation remains to be implemented 
for the lower-px DIS process (but is there in the higher-p^ direct ones). 

Therefore the current model should not be considered as the end of the road, but 
rather as providing a basic framework that could be further refined. In the end, data will 
have to tell whether the approach as such is viable or not. In the latter case, a simpler 
scenario would then be preferable, but more likely an even more complex one would be 
required. 
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